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^D [ Abstract 

We presenlp results of numerical modeling made for the galactic stellar 
and stellar-gas disk embedded in the spherical halo and bulge. The stellar 
disk is simulated by N-body system, the equations of hydrodynamics are 
1^ , solved by TVD-method. We used TREEcode-algorithm for calculation of 

^sj ' a self-gravity in stellar and gaseous components. The possibility of bars 

birth in a hot stellar disk because of gravitational instability of a cold 

^^ . gas component is investigated. The conditions of occurrence lopsided- 

\^ • galaxies from a axisymmetric disk as a result of gravitational instability 

r) I are explored. The self-consistent models of double bars are constructed 

• ■ and the dynamical stability of these structures is discussed. 
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1 Introduction 



Dynamics of many structures in disk galaxies is considerably determined by 
spherical subsystem properties and, in particular, by characteristics of density 
^ i distribution of dark halo in stellar disk limits. 

^SJ I The bar formation because of global bar-mode instability is impeded, if halo 

C^ . mass inside of optical radius surpasses disk mass M^ < (1 — 1.5) • Md [U [TO] . 

On the other hand, the observations data and N-body simulation give estimates 

Mh/Md > 1.5 for some galaxies [H [5]. Gas component is cold because of 

C"^ I radiative cooling and can be gravitationally unstable. The unstable modes in 

massive gas disk are capable to generate the bar even in the hot stellar disk in 
case of a high dispersion of stars velocities and at presence of the massive halo. 
A asymmetrical spiral structure (one-arm) and bar displacement concerning 
disk centre are typical distinctive features of a series SBcd-SBm galaxies (Mag- 
ellanic type). These properties are observed at LMC, NGC 55, 925, 1313, 1744, 



H ; 4490, 4618, 4625 etc. il US]. The formation mechanism of the displaced bar 

and other features of lopsided-galaxies can be caused by preferred growth of 
one-arms modes in gravitationally unstable disk and by subsequent interaction 
of these perturbations with a bar-mode at a nonlinear stage [12] . The late type 
galaxies contain more gas, than early type objects. Therefore question on influ- 
ence of gas on the bar displacement and asymmetry in disk structure requires 
of special study. 



1 2006ASSL..337..337K : Astrophysics and Space Science Library 2006, v.337, 
p. 337-344 // Astrophysical Disks / Alexei M. Fridman, et al (eds.). Publ. by 
Springer, Dordrecht, The Netherlands, 2006 
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Figure 1: Time dependences of the dispersion of stars radial velocities Cr (a) 
and sound speed of gas Cg (b) on various radiuses (curve 1 — disk center, 2 — 
disk periphery). 

The small-scale asymmetrical structures at disks center are very important 
for understanding of a phenomenon of nuclear galactic activity. The double 
bars can deliver gas to a active nuclei [TT]. The photometric data are the basic 
evidence about presence of the second inner bar at approximately 70 galaxies [S] . 
The self-consistent models of double bars were studied by the N-body method 
[21 [3] . Key problem of double bars is the question on dynamic stability of these 
systems. 



2 Modelling 

2.1 The numerical model of stellar-gas disk 

The 3D stellar disk simulation is based on A^-body model, taking into account 
an external field of the rigid matter distributions in bulge and halo. The gas 
disk model is constructed on the non- viscous equations of gasdynamics, and 
is complemented by gravitational forces on the part of stellar disk, spherical 
subsystem, and gas self-gravity also. 

The gas galactic disks are cold, as the sound velocity Cg is much less than 
the dispersion of radial velocities in stellar disks Cr- The radiative losses in the 
equation on energy are defined by quantity Q~: 



Q- ^A,_ 



{cl 



(c 
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■ Q^ for Cs > Csi [Q =0 for Cs < Csi) 



(1) 



Q — density, parameters A^.^ Csi, a, Cs2, l3 are free. The cooling of gas strongly 
grows in the case Cs — >■ Cs2, therefore restriction Cs < Cs2 is carried out always. 
We solved hydrodynamical equations by the method TVD-E. The self-gravity 
account in gas and stellar disks is based on TREEcode. We simulated disks 
with an exponential profile of surface density and radial scale L. 
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Figure 2: The gravitational instability of gas is the reason of bar formation in 
central area of the hot stellar disk. The isolines of a surface density are shown. 

If the stellar disk is on threshold of gravitational stability, small mass of 
cold gas component does not give to an additional heating of the stellar system 
(fig. [Ij. The gas mass is equal Mg — 0.08 • M^ inside the stellar disk in this 
model. And the massive halo {Mh = SMd) forbids the bar formation, as in 
stellar component (which besides is hot), and in gas disk. 

2.2 The stellar bar formation because of gravitational in- 
stability in the gas disk 

Let's consider models with halo mass Mh = (1-^2.5)- Md- In all cases the initial 
dispersion of star velocities and the massive halo provide gravitational stability 
of the stellar disk in absence of gas. The account of gas can qualitatively change 
evolution of system. 

The radiative cooling provides cold, gravitationally unstable state of gas 
component, it gives in formation in gas of non-axisymmetric structures, which 
in turn generate disturbances in the stellar disk (fig. [2|). There is the prompt 
bar formation in the stellar disk because of gravitational gas instability, if the 
models contain a lot of gas. It is important, that the instability of gas can 
generate the bar in hot stellar disk (fig. |3]), when Toomre's parameter exceeds 
Qt — Cr/cT > 2 in the region r < L. Such disks are stable without gas. 



2.3 Lopsided-galaxies 

Let's consider key influence of the gas component on effects of bars displace- 
ment and occurrence of asymmetry in isolated disk structure in the whole. The 
initial distribution of dispersion of stars velocities and the parameters of spher- 
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Figure 3: The isolines of surface density of stellar disk (left panel) and gaseous 
disk (right panel) after 5 rotation periods of outer part of disk. The standings 
of shock waves are well visible in region of the bar. 



ical subsystem suppose the slow formation of bar without the account of gas 
component. However the amplitude of one-arm harmonic {m = 1) is very small 
and formation of the lopsided-disk does not occur. 

The formation of harmonic to = 1 is possible in a massive cold gas sub- 
system. The nonlinear interaction of a one-arm mode and bar-mode (m = 2) 
in the gas disk is the reason of asymmetry of stellar disk also. The results of 
simulation in case of Mgas ~ 0.47 ■ Md in limits of r < 4L are shown in fig. 21 
With growth of relative gas mass we have amplifications of the bar displace- 
ment concerning centre of disk and power of spiral structure asymmetry. The 
considered mechanism of formation of lopsided-galaxies is most effective in case 
of small halo mass and if the halo scale exceeds the exponential disk scale in 2 
times and more. 



2.4 The problem of double-bars formation 

The bar formation requires not a hot initial disk and the halo mass Mhaio + 
Mbuige H'^Mdisk at halo scale a (1 — 4) • L (L — exponential disk scale). The 
birth of inner bar in numerical models occurs at presence of enough massive 
bulge {Mbuige ^0.3Md). In FiglSJ we show the distributions of surface den- 
sity logarithm lg{a) at the different time moments in model with M^ = M^, 
Mbuige = 0.6Mci, in which at particular stages there are structures such as 
double bars. 

The features of a kinematics of disk central region (r < 2L) at a stage of 
double bar are shown in a figEl The field of velocities in the stellar disk is 
the important information on existence of inner bar. The radial velocity U 
demonstrates four-areal structure, both for inner bar, and for primary bar. 
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Figure 4: The isolines of surface density of gas (a) and stellar (6) components 
are shown in the case Mt/Md = 1, Mgas/Md = 0.47. The initial Toomre 
parameter is equal Qt > 1-3 for the stellar disk in region r < 2L. 



The lifetime of double-bar does not exceed 1-2 rotation periods in the most 
ideal model of the stellar disk. The bending instability of disk and/or bar 
is the important factor of double-bars decay. The additional account of gas 
qualitatively changes result. The double-bars are not forming in self-consistent 
numerical stellar-gas models because of additional nonlinear perturbations. Ex- 
ternal asymmetrical potential (the tidal influence from the massive companion) 
gives similar result and such models do not give double-bars also. The con- 
clusion about a dynamically fast phase of existence of "double bar" is agreed 
with work [7], that the secondary bar not is real dynamically allocated struc- 
ture at observed galaxies, and represent a combination of objects with various 
morphology. 

3 Conclusions 



1. The bar formation in the hot gravitationally stable stellar disk can be gener- 
ated by the unstable cold gas disk. This mechanism generates the bar even in 
case of the massive halo Mh/Md ~ 1 — 2. 

2. The account of gas component strengthens the formation of asymmetrical 
structures (lopsided, mode to = 1) in a isolated disk as a result of gravitational 
instability in case of halo with small mass and large scale in comparison with a 
disk scale. At Magellanic type galaxies the relation of halo mass to disk mass 
in limits of optical radius on the average is less, than at systems of early types. 

3. The self-consistent models with the double bar are extremely unstable in 
relation to the various factors (transient spiral waves in a disk plane, bar warps, 
bending instabilities of a disk, tidal influence, gas component) at initial stages 
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Figure 5: The surface density distribution in different time (^, 5, 7, 8, 9 — 
double bars). 

of evolution. The conclusion about a very short phase of existence of systems 
such as "double bars" (transient nature) is made and similar structures can arise 
under special conditions at an initial stage of bar- mode development. 
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Figure 6: The distributions of the dispersions of velocity components (cr, Ctp, 
Cz) and the isohnes of azimuthal velocity V and radial velocity U . The areas 
of positive and negative radial velocity are shown by different shading for the 
main and inner parts of the stellar disk. 
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